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ABSTRACT 

We present the observational results of Galactic HII region Sli2-294, using 
optical photometry, narrow-band imaging and radio continuum mapping at 1280 

MHz, together with archival data from 2MASS, MSX and IRAS surveys. The 
stellar surface density profile based on 2MASS data indicates that the radius of 
the cluster associated with the Sh2-294 region is ~ 2.3'. We found an anomalous 
reddening law for the dust inside the cluster region and the ratio of total-to- 
selective extinction [Ry) is found to be 3.8 ± 0.1. We estimate the minimum 
reddening E{B — V) — 1.35 mag and distance of 4.8 ± 0.2 kpc to the region from 
optical color-color and color-magnitude diagrams. We identified the ionizing 
source of the HII region, and spectral type estimates based on different methods 
are consistent with a star of spectral type ~ BO V. The 2MASS JHKg images 
reveal a partially embedded cluster associated with the ionizing source along with 
a small cluster towards the eastern border of Sh2-294. The ionization front seen 
along the direction of small cluster in radio continuum and Ha images, might 
be due to the interaction of ionizing sources with the nearby molecular cloud. 
We found an arc shaped diffuse molecular hydrogen emission at 2.12 /im and a 
half ring of MSX dust emission which surrounds the ionized gas in the direction 
of the ionization front. The HIRES processed IRAS maps show two clumpy 
structures of high optical depth (tioo) and low color temperature (7(60/100)) 
at the eastern border of the nebula. Self consistent radiative transfer model of 
mid- to far-infrared continuum emission detected near small cluster is in good 
agreement with the observed spectral energy distribution of a B1.5 ZAMS star. 
The morphological correlation between the ionised and molecular gas, along with 
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probable time scale involved between the ionising star, evolution of HII region and 
small cluster, indicates that the star-formation activity observed at the border is 
probably triggered by the expansion of HII region. 

Subject headings: HII regions - star clusters - extinction - interstellar mat- 
ter - star formation process - infrared: ISM - ISM: individual (S 294) - radio 
continuum: ISM. 



1. Introduction 

Massive OB star formation takes place in the dense core of molecular cloud and are 
associated with dense stellar clusters. Their pre-main-sequence (PMS) life scales are much 
shorter compared to low mass stars. The stellar winds and radiation emitted from these 
stars eventually disrupt the molecular cloud, dissociate and ionize the gas and form the HII 
regions. The radiation from such stars can modify the interstellar grains present in their 
immediate vicinity, hence affect the interstellar extinction law (cf. Chini Sz Wargau 1994; 
Pandey et al. 2000). Thus, change in the interstellar extinction law affects the conclusion, 
such as photometric distances and age determination of the clusters from comparison with 
theoretical isochrones. All these short-comings therefore make it difficult to understand the 
formation and early phase evolution of such objects. The interaction of large radiation fields 
of massive stars with dense molecular gas produce photodissociation regions (PDRs), which 
is an important tool to understand the star formation, as evolution of nearby dense cloud is 
primarily governed by these interactions. The interaction of shock generated by the massive 
stars may also be responsible for triggering further star formation within the molecular cloud 
in which they are born (Elmegreen & Lada 1977). 

In this paper, we have studied Sh2-294/RCW3 region (ctaooo = 07''16'"34.5^ t^aooo = 
— 09°26'38") which was identified as an HII region (Sharpless 1959) and is part of Monoceros 
cloud complex. In the optical, it appears like a butterfly with two wings separated by dark- 
lane, roughly east-west. In the center of dark lane, a BO. 5 star is detected (Moffat, Jackson, 
& Fitzgerald 1979). Non-detection of H2O maser in the search by Henkel et al. (1986) 
for the Sh2-294 region indicates that the region is an evolved HII region. The radio map 
obtained by Fclli et al. (1978) at 4.99 GHz, shows good morphological agreement between 
the nebulosities in optical and radio upto a linear size of 7'. The molecular cloud associated 
with Sh2-294 region has been observed in CO (Blair, Peters, & Vanden Bout 1975; Blitz, 
Fich, & Stark 1982). The distance to Sh2-294 in the literature ranges from 2.5 to 4.6 kpc 
as quoted by FeUi & Harten (1981). The Sh2-294 region is coincident with an IRAS Point 
Source Catalog (PSC) source (IRAS 07141-0920) with an increasing spectrum from 12 to 
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100 yum. Wouterloot & Brand (1989) reported kinematic distance to the IRAS 07141-0920 
as 3.01 kpc and a luminosity of 3.6 x 10^ Lq for this IRAS source. 

To continue our multiwavelength investigations of massive star forming regions (cf. Ojha 
et al. 2004a; Ojha et al. 2004b; Ojha et al. 2004c; Tej et al. 2006), we have carried out 
a detailed study of Sh2-294 region as a useful example for understanding the massive star 
formation with the following layout of the paper. In Sect. 2, we describe our observations 
and the reduction procedures. In Sect. 3, we discuss other available datasets used in the 
present study. Sect. 4 discusses the stellar component and deals with the basic parameter 
determinations of the cluster associated to the massive star. Sect. 5 describes the interstellar 
matter (ISM) component and the physical environment associated to the massive star. Sect. 
6 is devoted for the discussion on star formation scenario in Sh2-294 region and in Sect. 7 
we present a summary of our results. 

2. Observations and Data Reduction 

2.1. Optical photometry 

UBV{RI)c CCD photometric observations were performed for the Sh2-294 region on 
2005 December 26 and 2006 March 06, using the 2K x 2K CCD system at the f/13 Cassegrain 
focus of the 104-cm Sampurnanand telescope (ST) of ARIES, Nainital (India). The 0.37 
arcsec pixel"^ plate scale gives a field of view (FOV) of 12'. 6 x 12'. 6 on the sky. The read-out 
noise for the system was 5.3 e~, while the gain was 10 e~ per ADU. To improve the signal- 
to-noise (S/N) ratio, observations were made in 2 x 2 pixel binning mode. Several bias 
frames and twilight flat field exposures were taken during the observation. We observed the 
standard area SAlOl (Landolt 1992) on 2006 March 06 several times during the night for the 
purpose of the determination of atmospheric extinction coefficients and for the photometric 
calibration of the CCD systems. The log of observations is given in Table 1. The CCD images 
were flattened using twilight sky fiats and bias subtraction was done using the IRAFj^ data 
reduction package. Then, for a given filter, frames of same exposure time were combined into 
one, to improve the S/N enabling us to detect fainter stars. The photometric measurements 
of the stars were performed using DAOPHOT II profile fitting package of MIDAS. The stellar 
images were well sampled and a variable PSF was applied from several uncontaminated stars 
present in the frames. The photometric errors at brighter magnitude level (V < 15) are 
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~ 0.01 mag, which increase with V magnitude and become large (~ 0.1 mag) at V ~ 22.5 
mag. The image parameters (e.g. chi & sharpness) and errors provided by DAOPHOT were 
used to reject poor measurements and only the stars having error in magnitudes < 0.1 are 
used in the further analyses. 

2.2. Optical spectroscopy 

To study further, a spectroscopic observation was made on 2006 March 16 for one of 
the sources, possibly responsible for the ionization of Sh2-294 region (cf. Moffat, Jackson, 
& Fitzgerald 1979 and sec §4.3). We used the Himalaya Faint Object Spectrograph Camera 
(HFOSC) available with the 2.01 m Himalayan Chandra Telescope (HCT), operated by the 
Indian Institute of Astrophysics, Bangalore (India). The instrument is equipped with a 
SlTe 2K X 4K pixel CCD. The central 2K x 2K region, used for imaging, corresponds to 
a FOV of ~ 10' X 10', with a pixel size of 0.296 arcsec. The sht spectra in the wavelength 
range, 3700-6800 A, were obtained using low resolution grism ^^7, with a slit width of ~ 
2". The spectroscopic observations were carried out under good photometric sky conditions. 
Bias and flat frames were obtained at the beginning and end of the observations. The one- 
dimensional spectra were extracted from the bias-subtracted and fiat-field corrected images 
using the optimal extraction method. Wavelength calibration of the spectra was done using 
FeAr lamp source. 

2.3. Optical and NIR Ncirrow-band imaging 

Optical CCD narrow-band images of the nebula around Sh2-294 region were obtained 
in Ha filter (A = 6565 A, AA = 80 A) and its nearby continuum filter (A= 6550 A, AA = 
80 A) using 104-cm telescope at ARIES, on 2006 September 28. Narrow-band near-infrared 
(NIR) observations were carried out in the rotational-vibrational line of molecular hydrogen 
- H2(1-0)S1 (A = 2.122 i^m, AA = 0.019 fim) and K continuum ( A = 2.264 /xm, AA = 0.054 
urn) with HCT on 2005 November 28 using the NIR imager (NIRCAM), which is a 0.8 - 2.5 
//m camera with 512 x 512 HgCdTe array. For our observations, the NIRCAM was used 
in the Camera-B mode which has a FOV of ~ 3'. 6 x 3'. 6. We obtained several dithered 
(by 20") frames of the target in order to remove bad pixels, cosmic rays and to eliminate 
the presence of objects with extended emission for construction of sky images. All images 
were bias-subtracted and fiat-field corrected in the standard manner using IRAF package. 
The final images are produced by subtraction of continuum images after aligning and PSF 
matching. 
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2.4. 



Radio continuum observation 



In order to probe the ionized gas component, radio continuum interferometric observa- 
tion at 1280 MHz was conducted on 2006 December 8, using tlie Giant Metrewave Radio 
Telescope (GMRT) array. Tlie GMRT has a "Y" -shaped hybrid configuration of 30 anten- 
nas, each 45 m in diameter. There are six antennas along each of the three arms (with arm 
length of ~14 km). These provide high angular resolution (longest baseline ~25 km). The 
rest of the twelve antennas are located in a random and compact 1x1 km^ arrangement 
near the centre and is sensitive to large scale diffuse emission (shortest baseline ~100 m). 
The observations near 1280 MHz were made using the full 16 MHz bandwidth. Details of 
the GMRT antennae and their configurations can be found in Swarup et al. (1991). The 
source, 3C147, was used as the primary fiux calibrator, while the source 0432-1-416 was used 
as a secondary calibrator to derive the phase and amplitude gains of the antennae. The data 
analysis was done using AIPS. The corrupted data (dead antennae, interference, spikes, 
etc.) were identified using tasks IBLED, UVPLT and VPLOT. The fiagging of the bad data 
was carried out using UVFLG and TVFLG. The image of the field was formed by Fourier 
inversion and cleaning algorithm (IMAGR). 



To complement our optical and NIR imaging data we have used the data for point 
sources around Sh2-294 region from the Two Micron All Sky Surve j0(2MASS) Point Source 
Catalog (Cutri et al. 2003) in the J (1.25 /im), H (1.65 /im) and K^-band (2.17 fim). The 
2MASS data have positional accuracy of better than 2". Our source selection was based on 
the 'read-flag' which gives the uncertainties in the magnitudes. To discuss the radial density 
profile of the cluster (see §4.1) associated with Sh2-294 region, we used 2MASS sources with 
'read flag' values of 1-6 from our sample. However to further improve the analysis, only 
sources with 'read-flag' values of 1-3 have been used to study the 2MASS color-color (CC) 
and color-magnitude (CM) diagrams (§4.6). It is also to be noted that we used 2MASS 
images for astrometric purposes of our optical images with matching tolerance of 1".0. 



■^This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project 
of the University of Massachusetts and the Infrared Processing and Analysis Center/ California Institute of 
Technology, funded by the NASA and the NSF. 
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Other Available Data Sets 



3.1. 



Near-Infrared Data from 2MASS 
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3.2. Mid- infrared data from MSX 

The Midcourse Space Experiment (MSX) surveyed the entire Galactic plane within 
|6| < 5° in four mid-infrared (MIR) bands centered at 8.28 /xm (A), 12.13 /xm (C), 14.65 
lim (D) and 21.34 fim (E) at a spatial resolution of 18". 3 (Price et al. 2001) and a global 
absolute astrometric accuracy of about 1".9 (Egan et al. 2003). The MIR data are useful 
to study the large scale distribution of warm (~ 1000 K) dust and emission from Polycyclic 
Aromatic Hydrocarbons (PAHs). The MSX A (50% peak response range: 6.8 - 10.8 fim) 
and C (11.1 - 13.2 /im) bands include several PAH features at 6.2, 7.7, 8.7, 11.3 and 12.7 
fim. Point sources have been selected from MSX Point Source Catalog Version 2.3 (Egan et 
al. 2003) and MSX A-band image is used to see PAH emission around Sh2-294 region. The 
flux densities of the MSX point source associated with Sh2-294 region are given in Table 2 
and used for constructing the spectral energy distribution (SED). 



3.3. Mid- and far-infrared data from /i?A5-HIRES 

The data from the Infrared Astronomical Satellite (IRAS) survey in the four bands 
(12, 25, 60 and 100 fim ) for the region around Sh2-294 were HIRES-processed (Aumann et 
al. 1990) at the Infrared Processing and Analysis Center (IPAC), Caltech, to obtain high 
angular resolution maps. Apart from determining the flux densities, the Ji?y4S'-HIRES maps 
have been used to generate the dust color temperature and optical depth maps of interstellar 
dust. The integrated flux densities of the IRAS PSC associated with Sh2-294 region are 
given in Table 2. 



4. Infrared and Optical Properties of Stellar Sources 

The Ha+continuum image and a false JHK^ color composite infrared image (J, blue; 
iJ, green; and K^, red ) of Sh2-294 region are shown in the left and right panels of Fig. 1, 
respectively. A close-up view of JHKg image reveals that there are two signiflcant amounts 
of stellar density enhancements in the region. The eastern region is situated at the border 
of Ho; emission, hence the region of interest of this study. 



■^This research made use of data products from the Midcourse Space Experiment. Processing of the data 
was funded by the Balhstic Missile Defense Organization with additional support from NASA Office of Space 
Science. This research also made use of the NASA/ IPAC Infrared Science Archive, which is operated by 
the Jet Propulsion Laboratory, Caltech, under contract with the NASA. 
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4.1. Radial Profile and Stellar Surface Number Density 

Sh2-294 is listed as a partially embedded infrared cluster in the catalog compiled by Bica, 
Dutra, & Barbuy (2003). To determine the boundaries of the cluster, we performed direct 
star counting, assuming spherical symmetry. We have selected sources which are detected in 
the 2MASS Kg band. To estimate the cluster radius we selected a region of 300" centered on 
IRAS 07141-0920 source. The center of the cluster was estimated by convolving a Gaussian 
kernel with the stellar distribution and taking the point of maximum density as the center. 
This was performed for both the axes to get the central coordinates of the cluster. The center 
was found to be at q;2ooo — 

07''16™33", ^2000 = -09°25'35". To determine the radial surface 
density, the Sh2-294 region was divided into a number of concentric rings from the cluster 
center. The projected radial stellar density in each concentric circle was obtained by dividing 
the number of stars in each annulus by its area. The error bars are derived assuming that 
the number of stars in each annulus follow Poisson statistics. To account for the contribution 
from the field stars, we selected a control field (0:2000 = 07'*17"*33^ ^2000 = -09°25'24") which 
is ~ 15' to the east of IRAS 07141-0920 source. Figure 2 shows the radial density profile of 
the cluster (open circles) and the horizontal dashed hue shows field star density estimated 
from the control field. A clear gradient in surface number density distribution of Sh2-294 
region confirms the existence of clustering around the center of the region. Several studies 
have shown that the radial density profiles of young embedded clusters can be fitted by both 
the inverse radius (r~^) model and the King's model (e.g. Horner et al. 1997; Teixeira et al. 
2004; Baba et al. 2004). Neglecting the tidal radius, the King's profile can be written as 



where is the core radius (the radius at which the surface density falls to half of the 
central density /o) and a is the constant for background offset due to field stars. 

For the Sh2-294 cluster, we found that the King's model has a better fit (reduced 
= 0.96) as compared to 'r~^' model (reduced = 1-04). Within errors, the King's profile 
(dotted curve in Fig. 2) merges with the background field at ~ 2.3', which corresponds 
to ~ 3.2 pc at a distance of 4.8 kpc (see §4.5). We take this as the cluster radius. The 
background level as estimated from the control field is ~ 4.7 stars arcmin"^ which is in 
reasonable agreement with the value of ~ 3.8 stars arcmin~^ obtained from King's profile. 

The contour map of the stellar surface number density (SSND) around Sh2-294 region 
was made using kernel method (Gomez et al. 1993; Silverman 1986) and is shown in Fig. 3. 
The stellar density distribution is smoothened by 12" x 12" samphng box. At the distance of 
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4.8 kpc, the unit cell corresponds to ~ 0.28 pc x 0.28 pc. The contour levels are drawn from 

3 times of the background level with an increment of 2 stars/pc^, where the background 
level is taken as 4 stars/arcmin^ (2 stars/pc^). The peak of the projected stellar density 
distribution was found to be 02000 = 07'^16™33'', ^2000 = — 09°25'29", which is close (6") to 
the position found by convolving a Gaussian kernel as discussed above. The King's profile 
fitting yields a central density of 44 ± 13 stars/arcmin^ (~ 22 stars/pc^), correlates well 
with ~ 22 stars/pc^ of the peak of SSND. The clustering in the central region is more 
obvious in the surface density map. However, there is evidence that Sh2-294 possesses 
distinct substructures. The substructures in the east and northwest directions indicate the 
density enhancement, which can also be seen as a second peak in radial density distribution 
around ~ 1.2' from the central peak (see Fig. 2). The eastern substructure is roughly 
spherical in shape and more prominent with peak approximately at CI2000 — 

07^16^38' and 

^2000 = — 09°25'40", having central density of 13 stars/pc^. This substructure is at a distance 
of ~ 1.23' from the cluster center of Sh2-294 region and is also located beyond the faintest 
nebulosity associated with the central cluster seen in the 2MASS Kg band image. Therefore, 
it might be an other independent cluster but still embedded in the molecular cloud complex 
harboring Sh2-294. Hereafter, we refer to this region as 'A' for our further discussion in the 
paper. In Figure 1 (right) the 'region A' is represented by a white circle in JHKg color 
composite image. 



4.2. Probable members of the cluster 



To study the various parameters of the cluster, it is necessary to remove field star 
contamination from the sample of stars in the cluster region. Membership determination 
is also crucial for assessing the presence of PMS stars, because both PMS stars and dwarf 
foreground stars occupy similar positions above the zero-age main sequence (ZAMS) in the 
color-magnitude (CM) diagram. In the absence of proper motion, we used statistical criterion 
to estimate the number of probable field stars in the cluster region by selecting a control 
region of same area from our CCD image (control field is used to estimate the field star 
contamination). Figures 4a and 46 shows V, {V — I) CM diagrams for the cluster and control 
field regions, respectively. The contamination due to background field population is clearly 
visible in the cluster region CM diagram. To remove contamination of field stars from the 
cluster region sample, we statistically subtracted the contribution of field stars from the CM 
diagram of the cluster region following the procedure given by Mateo Sz Hodge (1986). For 
a randomly selected star in the V, {V — I) CM diagram of the field region, the nearest star 
in the cluster's V, {V - I) CM diagram within V ± 0.25 and {V - I) ± 0.125 of the field star 
was removed. The statistically cleaned V, {V — I) CM diagram of the cluster region is shown 
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in Fig. 4c which clearly shows the presence of PMS stars in the cluster. The presence of a 
few field stars in the statistically cleaned CM diagram can still be noticed. These probable 
field stars are identified by using the CC diagram and marked as triangles in Fig. 4c as well 
as in Fig. 6 (see §4.4). 

4.3. Spectral type of ionizing source 

The brightest member found in our statistically cleaned CM diagram is situated at the 
center of optical nebula having 0:2000 = 07'*16'"33*, ^2000 = 

-09°25'22". The photometric 

spectral type of this star is of BO. 5, with V = 13.99 and B — V = 1.03, which is the source 
of ionization in Sh2-294 (Moffat, Jackson, & Fitzgerald 1979). This is in agreement with our 
measurements, V = 14.04, B — V = 1.04 and My = -3.56 at the assumed distance 4.8 kpc 
(see §4.5), correspond to BO. 5 V spectral type (Schmidth-Kaler 1982; Panagia 1973). The 
second brightest member of the cluster has V — 15.41 and My = -2.20, which is consistent 
with spectral type between B2-B3 V (Schmidth-Kaler 1982; Panagia 1973). As stars of 
spectral type B2 or earlier are known to produce a significant proportion of their radiation 
as photons with energy high enough to form HII region, therefore the star situated at the 
center is most probably the ionizing source of Sh2-294 HII region. 

To determine the spectral type of this ionizing source, we extracted ID-spectrum (see 
§2.2), with aperture centered on the star. Figure 5 shows the spectrum of the star in the 
range 3950-4750 A. Besides hydrogen Balmer lines (H^, H^, H^, etc), the spectrum displays 
Hell line at A ~ 4200 A, and weak Hell hue at 4686 A, along with OH/CHI blend at 
4650 A, indicates that the spectral type of the star is earlier than Bl. Specifically Hell 
line at A 4200 A is seen in stars of spectral type B0.5-B0.7 or earlier. On the other hand, 
the weakness of Hell line at A 4541 A relative to Hel line at 4471 A indicates an spectral 
type latter than 09, while Hell A4686 A < Hel A4713 A favors the classification of ~ BO. 
Comparing the composite spectrum with that in the atlas of low-resolution stellar spectra 
of Jacoby et al. (1984) and with the atlas of Walborn & Fitzpatrick (1990), we propose that 
the exciting source for Sh2-294 has an equivalent spectral type of BO ± 0.5 V. 

4.4. Extinction 

The extinction in star cluster arises due to two distinct sources; (i) the general ISM in 
the foreground of the cluster {E{B — V)jnin): and (ii) the localized cloud associated with 
the cluster {AE{B -V)^ E{B - V)^ - E{B - V)min), where E{B - V)^ is the reddening 
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of the star embedded in the parent cloud. The former component is characterized by the 
ratio of the total-to-selective extinction Ry (= Ay/E[B — V))= 3.1 (Wegner 1993; He et 
al. 1995; Winkler 1997), whereas for the intracluster regions of young clusters embedded 
in dust and gas cloud the value of Ry varies significantly (Tapia et al. 1991; Pandey et 
al. 2003). However, Hll regions associated with a large amount of gas and dust often show 
larger values of Ry (Forte 1978; The & Groot 1983; Chini & Wargau 1990; Pandey et al. 
2000). 

In the absence of spectroscopic observation, the interstellar extinction E{B — V) towards 
the cluster region can be estimated using the {U — B)/{B — V) CC diagram. The CC 
diagram of the cluster region is presented in Fig. 6, which contains data points for all the 
stars having V < 18 mag. The figure clearly shows a large amount of contamination due 
to foreground stars. The probable field stars in the statistically cleaned CM diagram follow 
the sequence of foreground stars and are marked as triangles in Fig. 6. The probable cluster 
members manifest a more reddened sequence than the foreground sequence. To match the 
observations, the ZAMS by Schmidt-Kaler (1982) is shifted along the reddening vector having 
a normal slope of = 0.72. The CC diagram indicates that the foreground population 

is reddened by E{B — V) = 0.2 mag. The minimum reddening for the cluster region is 
estimated to be E{B — V)min = 1-35 mag. 



4.4- 1- Ratio of Total-to-Selective Extinction 

The variable amount of absorbing grains in the line of sight determines the total extinc- 
tion towards each star. In practice, the color excess for a star is obtained by comparing its 
intrinsic color with the observed color, and the total visual extinction is derived by using the 
shape of dust-reddening law. The direct method, however, cannot be easily applied because 
the intrinsic colors of most of the stars are usually unknown, so that an indirect way of sta- 
tistical approach can be adopted on the basis of number of stars detected in the investigated 
region. It has been shown by Chini and Wargau (1990) that two color diagrams (TCDs) of 
the form of {V — A) vs. {B — V), where A is one of the wavelengths of the broad-band filters 
(R, I, J, H, K), provide an effective method of separating the abnormal extinction arising 
within such clusters from the normal extinction produced by general ISM. Chini & Wargau 
(1990) and Pandey et al. (2000) used TCDs to study the anomalous extinction law in the 
clusters M 16 and NGC 3603, respectively. 

To construct TCDs, we used the probable cluster members obtained from statistically 
cleaned CM diagram (cf section 4.2). The stars showing infrared excess are not included in 
the analysis. Figure 7 shows the {V — X)/{B — V) TCDs for all the probable members of the 



cluster region. The slopes of the distribution are obtained by using the weighted least-square 
fit to the data points. To derive the value of Rduster (ratio of total-to-selcctivc extinction in 
the cluster region) we use the approximate relation (cf. Neckel & Chini 1981), 



p _ "^cluster T3 (r^\ 

^duster ^ ^normal V J 

^normal 

Here mduster and rrinormai are the slopes of the distribution followed by cluster stars 
and normal MS stars (taken from the stellar models by Bertelli et al. 1994) respectively. 
The slopes are given in Table 3. The Rnormai represents to the value of 3.1. The weighted 
mean value of Rduster comes out to be 3.8 ± 0.1, indicating an anomalous reddening law for 
Sh2-294 region. Here we would like to mention that the value of Rduster obtained by using 
equation (2) should be considered as an approximate estimation of R in the cluster region. 

To further test the reddening and Ry value, we used the observed properties of ionizing 
source and assumed that the spectral type of the star is BO V. The observed V magnitude 
of the ionizing star is 14.04 with B — V — 1.04. Comparing with intrinsic colors of Johnson 
et al. (1966), the reddening E{B — V) to the star corresponds to 1.35 mag. The Ry value 
to the star can be calculated using the empirical relationship established by Whittet (1976) 
between Ry, the color excesses E{B — V) and E{V — K), 



Rv^l.lQ^E{V -K)/E{B-V), (3) 

To derive Ry for the star, the JHKs magnitudes of the star arc transformed from 2MASS 
system to Koornneef system using the relation given by Carpenter (2001) and intrinsic V — K 
color obtained from Koornneef (1983). The Ry value thus obtained is ~ 3.9. Hence, the 
reddening and ratio of total-to-selective extinction estimated as above agree quite well with 
the value found from TCDs. 



4.5. Distance and age of Sh2-294 

The ages of young clusters are typically derived from post-main-sequence evolutionary 
tracks for the earliest members if significant evolution has occurred or by fitting the low-mass 
contracting population with theoretical PMS isochrones. The fitting of theoretical zero-age 
main sequence (ZAMS) to the observed sequence yields distance. The distance of the cluster 
associated to Sh2-294 is derived using statistically cleaned CM diagrams as shown in Fig. 
8. Using the foreground reddening E{B — V) — 1.35 mag towards the cluster region and 
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E{V - R)/E{B -V) ^ 0.60, we visually fitted the ZAMS (sofid curve) given by Schmidt- 
Kaler (1982) and Walker (1985) to the V, B-V and V, V-R CM diagrams respectively. The 
visual fit yields a distance modulus of (m — M)v — 17.6 ± 0.10, which corresponds to a 
distance of 4.8 ± 0.2 kpc. Our distance estimation is in agreement with the most widely used 
distance of 4.6 kpc for this source. Theoretical isochrone for 4 Myr (Z=0.02) by Bertelli et 
al. (1994) is also shown in Fig. 8. Since massive stars do not show significant evolution, a 
maximum age of 4 Myr can be expected for the cluster. In Fig. 8 we have also overplotted 
PMS isochrones from Siess, Dufour, & Forestini (2000) for 1 and 5 Myr. The open circles 
with dots inside represent the optical counterpart of NIR excess stars (see §4.6.1), which are 
the probable PMS sources of the cluster. The positions of these sources in the CM diagram 
indicate an age spread of 1-5 Myr. 

4.6. Nature of sources associated to Sh2-294 

In JHKs color composite image (see Fig. 1 right), one could see that the 'region A' 
is significantly redder than other parts of the region. The sources in this region have no 
optical counterpart implying that they are still deeply embedded in the molecular cloud of 
high visual extinction. 

4-6.1. NIR Color- Color Diagram 

In order to examine the nature of the stellar population, we use (H — K) vs. (J — H) CC 
diagrams of the Sh2-294 star-forming region and the control field. In Fig. 9 we have plotted 
141 and 79 sources for the cluster region and the control region, respectively, having 2MASS 
'read-flag' values of 1-3 (good quality photometric magnitudes) in all the three JHKs bands. 
To plot the stars in the CC diagrams, we transformed the magnitudes and colors to the CIT 
system. The solid thin curve represents the locus of MS and thick-dashed curve is for giant 
stars taken from Bessell & Brett (1988). The dotted line represents the classical T Tauri 
(CTT) locus as determined by Meyer, Calvet, & Hillenbrand (1997). The thick continuous 
solid curve represents the Herbig Ae/Be locus (Lada & Adams 1992). Presuming universal 
nature of the extinction law in NIR region (Mathis 1990), we plotted reddening vectors of 
the normal reddening law as parallel dashed lines. We assumed that Aj/Ay — 0.265, Ah /Ay 
= 0.135, Ak/Av = 0.090 for CIT system (Cohen et al. 1981). 

The arrow in Fig. 9a points to the position corresponding to the central brightest star 
in the cluster. We calculated the spectral type of the star by shifting it along the reddening 
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vector to the intrinsic main sequence locus. We found the spectral type of the brightest star 
in the cluster to be ~ BO. In Figure 9a we have shown the sources belonging to 'region A' 
in squares; these are the sources in a circular area of radius 20" around the second peak 
(eastern substructure) seen in the SSND diagram (see Fig. 3). Figure 9a shows that the 
average visual extinction (Ay ~ 7.4 mag) for the sources of 'region A' is higher than the 
average Ay (~ 3.3 mag) of the sources associated to other parts of Sh2-294 region. We 
have calculated the extinction by de-reddening the stars of "F" region (see below) in the CC 
diagram to the M0-M6 part of main-sequence locus. The high Ay indicates that the stars 
in 'region A' are more deeply embedded thus suffers a higher extinction. 

We have classified the CC diagram into three regions (e.g Sugitani et al. 2002; Ojha 
et al. 2004a,b) to study the nature of sources. The "F" sources are located within the 
reddening band of MS and giant stars. These stars are generally considered to be either 
field stars, Class III objects or Class II objects with small NIR excess. "T" sources are 
located redward of region "F" but blueward of the reddening line projected from the red 
end of the CTT locus of Meyer, Calvet, & Hillenbrand (1997). These sources are considered 
to be mostly T Tauri stars (Class II objects) with large NIR excess. "P" sources are those 
located in the region redward of region "T" and are most likely Class I objects (protostcllar 
objects). There may be an overlap in NIR colors of upper-end band of Herbig Ae/Be stars 
and lower-end band of T Tauri stars in the "T" region (Hillenbrand et al. 1992). 

Majority of sources in our sample arc distributed in the "F" region. A total of 36 sources 
show NIR excess (i.e. sources right to the "F" region), out of which 18 sources arc populating 
in the "T" and the "P" regions, characteristic of young stellar objects (YSOs) of Class II and 
I. However, this is the lower limit as several cluster members detected in the Ks band, may 
not be detected in the other two shorter NIR wavelength bands. The fraction of NIR excess 
stars in a cluster is an age indicator because the disks/envelopes become optically thin with 
the age. For young embedded clusters with age < 1 x 10^ yr the fraction is ~ 50% (Lada 
et al. 2000; Haisch, Lada, & Lada 2000) and the value is ~ 40%, as in case of the YSOs in 
the Taurus dark clouds, which have an estimated age of about ~ 1 — 2 x 10^ yr (Kenyon & 
Hartmann 1995). The fraction decreases to ~ 20% for the clusters with age ~ 2 — 3 x 10^ yr 
(Haisch et al. 2001; Teixeira et al. 2004). We have detected 6 point sources in 'region A' in 
all the three bands within the 2MASS detection limits. One of the sources located left of the 
reddening vector of giant stars (leftmost reddening vector in Fig. 9a) is probably a field star 
of the region whereas two sources show NIR excess, characteristic of YSOs of Class II. The 
fraction of the NIR excess stars (2/5, 40%) in 'region A' thus suggests an upper age limit 
of ~ 1 — 2 X 10^ yr of the cluster. We also found three 2MASS sources detected only in Kg 
band within 'region A' (see §6.1), which could be YSOs in their earliest stages still embedded 
in the cloud. The non- detection of these sources in J and H bands limits our analysis to 
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examine the exact nature of each source. If we presume that these sources are YSOs having 
circumstellar disks, then the fraction of NIR excess stars (5/8, 62%) would suggest an age of 
< 1 X 10^ yr to 'region A'. However, to further improve the statistics, observations in L and 
M-band are needed. In the CC diagram of control field, most of the sources are spread over 
lower part of "F" region; they are mainly foreground stars. In comparison the distribution 
of sources in the cluster region shows a mixture of field stars and reddened sources probably 
having warm circumstellar dust, characteristic of young PMS objects (Lada & Adams 1992). 
In order to examine the statistical properties of the cluster sources, the histograms oi H — K 
colors are shown in Fig. 10. The color distribution of sources belonging to control field is 
shown in shaded area, whereas, for the cluster region the distribution is shown in thick lines. 
The histogram shows two peaks in the range 0.05 < H — K < 0.2 and 0.35 < H — K < 0.5 
for the cluster region. The peak oi H ~ K color distribution of control region is ahead of 
cluster region in the range 0.05 < H — K < 0.2. Hence, probably the first peak oi H — K 
color distribution of cluster region may be due to foreground population of the region. On 
the other hand, a large number of sources were found in the cluster region with H — K > 0.4. 
The existence of these reddened sources indicates a population of young stars still embedded 
in significant quantities of intervening gas and dust. 



4-6.2. NIR Color- Magnitude Diagram 

Further information about the nature of the sources can be extracted from the (J — H) 
vs. J and {H — K) vs. K CM diagrams. In [H — K) vs. K CM diagram (Fig. 11) we have 
plotted all the 141 sources detected in the JHR^ bands plus additional 11 sources which are 
detected in the H and Ks bands only. The vertical solid lines represent ZAMS curves (for a 
distance of 4.8 kpc) reddened by Ay = 0, 10 and 20 mag, respectively. The parallel slanting 
lines represent the reddening vectors to the corresponding spectral type. An apparent MS 
track is noticeable with Ay ~ 1.5 mag in this diagram, however a comparison of it with a 
similar diagram for the stars in the control field shows that it is a false sequence caused by 
field stars in the foreground of Sh2-294 region. The foreground sources in cluster region are 
oi H — K < 0.2. YSOs (Class II and I) found from the CC diagram (Fig 9a) are shown as 
star symbols and filled triangles, respectively. The probable Herbig Ae/Be sources are shown 
as crosses and the sources belonging to 'region A' are marked in squares. Sources detected 
only in H and Kg bands are represented by large filled circles which may also be YSOs. For 
the sources detected in H and Kg bands, we have taken the Kg magnitude as Kcit due to 
absence of appropriate transformation equation between K2MASS to Kcit system but the 
colors are transformed to the CIT system. Out of 11 sources detected in H and Kg bands, 
only 6 of them show H — K > 1.0, hence they are probably young stars in their earliest 
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evolutionary phases. 

The mass of the probable YSO candidates can be estimated by comparing their locations 
on the CM diagram with evolutionary models of PMS stars. We used {J — H) vs. J CM 
diagram instead of the {H — K) vs. K color-magnitude diagram to minimize the effect of 
the "excess" emission in the NIR to derive the ages and masses of the YSOs. Figure 12 
represents intrinsic (J — H) vs. J CM diagram for probable YSOs (Class II and I). The 
symbols are same as in Fig. 11. To produce the intrinsic CM diagram, the extinction in 
front of YSOs are derived by tracing back along the reddening vector to the CTT locus or 
its extension in [H — K) vs. (J — H) CC diagram (Fig. 9a). Here we would like to mention 
that the present derivation of extinction is an approximate approach. The thin solid line 
represents PMS isochrone of 1x10^ yr from Siess, Dufour, & Forestini (2000), at an assumed 
distance of 4.8 kpc. Arrow hues are the reddening vectors for sources with masses of 0.5 
and 3.5 M© respectively for 1x10^ yr. As can be seen in Fig. 12 the majority of YSOs have 
masses in the range of 0.5 to 3.5 M© indicating that these may be T-Tauri stars. 



5. Description of the HII region 

5.1. Emission from Ionized Gas 

Figure 1 shows the Ha+continuum image of Sh2-294 region. The nebula looks loose 
bipolar in nature, having an elongated obscuring dark-lane extending approximately in the 
east-west direction cutting across at the center. The obscuring dark lane at the center 
suggests the presence of a molecular cloud. Indeed Blair, Peters, & Vanden Bout (1975) 
detected CO (J = 1 — >• 0) emission over a 6' x 6' area with peak in CO intensity at q;2ooo = 
07'^16'"32^ ^2000 = -09°25'23". The hues were centered at Ylsr = +32 km s'^. The dark- 
lane divides the nebula into two parts, with bright nebulosity to the south and less bright to 
the north, which decreases in surface brightness away from the obscured region. The southern 
component is brighter and has sharp arc structure in the east direction giving an impression 
of possible ionization front (IF). Whereas the northern component is more extended and 
very filamentary in structure, it gives a stunning view of the dust content of the region, 
proving that a significant amount of dust is concentrated in the string-like filaments across 
the region. This structured feature suggests an inhomogeneous spatial distribution of the 
interstellar matter. 

Felli & Churchwell (1972) reported the radio continuum emission from Sh2-294 region in 
their low resolution (HPBW = 10') 1.4 GHz survey of Galactic HII regions. Because of this 
resolution the object appeared featureless. However, the observation with VLA (Fich 1993) 
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with resolution of 40" shows that there is a gradient in free-free emission surface brightness 
from southeast to northwest part of the nebula. Figure 13 represents the radio continuum 
image of Sh2-294 with GMRT at 1280 MHz. The radio continuum image has a resolution of 
33" X 17". The map has an rms noise of ~ 0.26 mJy/beam. The integrated flux density in 
the map is of ~ 0.4 Jy, obtained by integrating down to 3(T level, where o is the rms noise 
in the map. The radio morphology shows that the thin ionization-bounded zone in south- 
eastern direction is an indication of IF. The decreasing intensity distribution on the opposite 
side signifies the region could be density bounded in northwest direction. From the optical 
and radio morphologies, it appears to be the flow of ionized gas in northwest direction, 
as expected in case of champagne model (Tenorio-Tagel 1979). Figure 14 represents the 
radio contours of GMRT at 1280 MHz superimposed on the band image of 2MASS in 
logarithmic scale. The bright source which is the possible source of ionization is located 
at the geometric center of the radio map and the embedded sources of the 'region A' are 
seen in the direction of IF, which could be the probable young PMS sources of the Sh2-294 
region. The radio emission peaks at 0:2000 = 07''16'"37*', ^2000 = — 09°25'45". At the peak of 
radio emission, we did not see any significant amount of Ha emission, which could be due 
to presence of dense ISM/cloud which obscures the optical emission. 

Assuming that Sh2-294 is roughly spherically symmetric ionization bounded, and, ne- 
glecting absorption of ultraviolet radiation by dust inside the HIT region, the measured flux 
density of 0.4 Jy, together with assumed distance, allow us to estimate the number of Lyman 
continuum photons {Nhyc) emitted per second, and hence the spectral type of the exciting 
star. Using the relation derived by Lequeux (1980) and assuming an electron temperature of 
8000 K, we estimate logA^L^c = 47.90, which corresponds to a MS spectral type earlier than 
BO (Panagia 1973) and close to BO. 5 (Vacca, Garmany, & Shull 1996; Schaerer & de Koter 
1997). respectively. The differences in the spectral type is expected due to uncertainties in 
the effective temperature and ionizing fluxes of massive stars, for a given spectral type. 



5.2. Dust continuum emission 

The strong mid- to far-infrared continuum emission detected from the local region of the 
interstellar cloud near 'region A' has been modelled using a radiative transfer scheme (Mook- 
erjea & Ghosh 1999). The energy from the central embedded source is transported through 
the dust and gas (only hydrogen) components of the spherical cloud and the emergent spec- 
trum is compared with the observed SED (see Table 2). The best fitting model is compared 
with the IRAS and MSX data in Fig. 15 (solid curve). The far infrared emission has been 
integrated over a 3' (diameter) region centered on IRAS 07141-0920 from the 7i?A,S-HIRES 
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maps. There is however an excess emission at 12 //m, which could be due to a component of 
very small grains or polycyclic aromatic hydrocarbons (Degioia-Eastwood 1992). The fit is 
reasonably good. The parameters of this model are : r° density law, radial optical depth of 
5.0x10"^ at 100 fim, central dust-free cavity size of 0.012 pc, total luminosity corresponding 
to a ZAMS star of type B1.5. It may be noted that the predicted radio emission at 1280 
MHz is too low compared to that observed locally (170 mJy within 3' diameter). With the 
aim of reproducing the observed radio emission, another model has been constructed. In 
order to fit the radio data well, an entirely different type of the embedded ZAMS star (06.5) 
was necessary. The resulting emergent spectrum in this case is also presented for compar- 
ison in Fig. 15 (dashed curve) which is at variance with the observed infrared emission. 
Hence, wc propose that while the continuum emission from the dust component is powered 
by the embedded source, which might be so heavily obscured that it emits mainly in mid- 
to far-infrared, the ionization of the gas is maintained externally (the BO star discussed in 
§5.1). 



5.3. Photodissociation region 

Figure 16 shows the continuum-subtracted B.2{i') = 1 — 05(1) 2.12 fim image smoothed 
with a 3 X 3 boxcar. The image has a field of view of 3.6' x 3.6'. Due to low signal-to-noise 
ratio, we can roughly describe only some of the characteristics of H2 emission. There is 
significant amount of diffuse H2 emission seen at the center of the nebula, where the dark- 
lane is located, which is correlated well with the CO peak. This kind of structure has been 
observed in the star-forming region S255-2, where the central core of S255-2 is immersed in 
a region of H2 emission with marginal evidence for limb brightening to the south (Howard, 
Pipher & Forrest 1997). We also found an arc-shaped diffuse H2 emission feature towards 
cast, as seen in case of Sh2-219 (Deharveng et al. 2006) in the direction of IF, which could 
be due to the interaction of massive star(s) with the nearby molecular cloud. The diffuse 
H2 emission can be excited by Lyman and Werner UV photons (e.g., Chrysostomou et al. 
1992) from the massive stars and thus can trace photodissociation (PDR) regions; but for a 
dense gas {hh ~ 10^ cm~^), H2 emission alone cannot be used to identify PDR and hence 
another PDR tracer hke Polycyclic Aromatic Hydrocarbon (PAH) emission is necessary to 
trace PDR around nebula (Giard et al. 1994). 

The far-UV radiation (6 < hu < 13.6 eV) can escape from HII region and penetrate 
the surface of molecular clouds, leading to the formation of PDR in the surface layer of the 
clouds. PAHs within the PDR are excited by the UV photons re-emitting their energy at MIR 
wavelengths, particularly between 6-10 iim. The A-band (8 iim) of MSX includes several 
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discrete PAH emission features (e.g. 6.6, 7.7 and 8.6 /xm) in addition to the contribution 
from thermal continuum component from hot dust. Figure 17 shows Ha image superimposed 
with MSX A-band contours. The 8 /xm emission displays partial ring in eastern side of Sh2- 
294 region, with diffuse emission extending towards the dark lane of optical nebula. The 
ring of PAH emission lies beyond the IF, indicating the interface between the ionized and 
molecular gas. The absence of 8 yum emission in the interior of HII region is interpreted as 
the destruction of PAH molecules by intense UV radiation of the ionizing star. The peak of 
PAH emission matches closely with diffuse Hg emission seen in the east. While comparing 
the PAH emission with Ha emission, it appears that the diffuse emission of MSX corresponds 
to dark optical features of Ha emission. 



5.4. Dust optical depth and temperature distribution 

The temperature of interstellar dust is determined by a balance between radiative heat- 
ing and cooling processes. Heating occurs primarily by absorption of stellar UV photons, 
and cooling by re-radiation in mid- and far-infrared. A dust color temperature of each pixel 
using the flux ratio at 60 and 100 /im can be determined by assuming a dust emissivity law. 
We have used HIRES processed IRAS maps to generate maps of the dust color temperature 
(T(60/100)) and optical depth (r25 and rioo) around the Sh2-294 region. The intensity maps 
at 12, 25, 60 and 100 fim were spatially averaged before computing T(60/100) and (r25 ,rioo) 
in a similar manner as described by Ghosh et al. (1993) for an emissivity law of eA oc A~^. 

The HIRES processed IRAS maps at 12, 25, 60 and 100 /xm are shown in Figure 18. The 
infrared emission in 60 and 100 /xm maps peaks at ^2000 = 07^16™'37*, ^2000 = — 09°25'51". 
The peak of infrared emission coincides (~ 9") with the 1280 MHz GMRT radio continuum 
peak and lies ~ 16" away from the peak of central density found in SSND distribution 
(see §4.1) for 'region A'. The maps at 60 and 100 /xm trace the distribution of colder dust. 
The peak of cold dust, along with the peak of radio continuum suggest that the gas-dust 
coupling is rather efficient in the 'region A', which is normally the case for regions of recent 
star formation. The dust emission seems to favor southeast to northwest orientation as for 
the cloud of Sh2-294 region. The extensions of the 60 and 100 dust emission maps have 
diameters larger than 12 and 25 fim maps. This is due to the fact that the former has its 
main contribution from the cooler component of dust grains lying in the outer envelope of 
the cloud. Figure 19 shows the contour maps of dust optical depth at 25 & 100 fim. The 
Tioo map represents low optical depth at the position of MSX point source and/or 'region 
A', however two peaks (with a peak value of tioo ~ 0.001) can be seen at the positions 
along northeast and southwest of 'region A', indicating presence of high densities in the 
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eastern side of Sh2-294. A gradient is seen in the optical depth map from east to northwest 
direction. The low optical depth in the northwest direction may be due to the lower dust 
column density along the path. Figure 20 presents a contour map of T(60/100). The dust 
color temperature varies from 30 to 40 K in the map. The T(60/100) and rioo distributions 
suggest an anticorrelation, such as the peak position of rioo map corresponds to a low color 
temperature of ~ 31 K. The relatively high temperature around 'region A' may probably be 
due to the radiation from young stars of the region. The cometary morphology of the map 
could be due to inhomogeneous distribution of dust in Sh2-294 region. 



6. Star formation scenEirio eiround Sh2-294 
6.1. Spatial distribution of young sources 



Figure 21 shows the continuum-subtracted molecular hydrogen 2.12 //m image of field 
of view 3.6' x 3.6', overlaid by the GMRT radio continuum contours (blue color). The 
red contours represent the MSX 8 fim emission, which arc generally attributed to PAHs. 
The spatial distribution of Class II, Class I, probable Herbig Ae/Be sources and sources 
detected only in H and Kg bands are shown in different symbols. The symbols for the 
sources are same as in {H — K) vs. K CM diagram (Fig. 11). We also find that there 
are five sources detected only in the band but not in H and J bands. These sources 
may also be probable YSO candidates; we represented them by red filled circles in Fig. 
21. We found a non-uniform distribution of young sources in Sh2-294 star-forming region. 
The complex seems to be composed of sources of different evolutionary stages aligned from 
northwest to southeast along the dark lane seen in Hct image (see Fig. 1). They are mainly 
concentrated around the region where diffuse H2 emission is seen. Most of the Herbig Ac/Be 
and Class II sources are seen in northwestern direction but in contrast, sources detected in 
H and K^^ bands and only in K^^ band are seen towards southeastern direction. The spatial 
distribution suggests that the southeastern region is less evolved than other parts of Sh2-294 
star forming region. The regions (northeast and southwest) which are most probably ionized 
by a massive star located at the center of the nebula, have lack of Class II, Class I and other 
probable YSOs, which implies that star formation activity has subsided in these directions 
due to intense UV photons. We identified a source from MSX PSC at q;2ooo = 07'*16"*38*, 
(^2000 = — 09°25'40" which coincides with the peak by found SSND distribution for 'region A'. 
In Sh2-294 region the young sources in and around 'region A' are observed in PDR region 
along the direction of IF. The morphology of the PDR and II2 emission in comparison with 
radio continuum emission gives evidence of interaction of massive star(s) with nearby cloud. 
The close association of the partially embedded cluster ('region A') adjacent to HII region in 
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the direction of IF therefore suggests that the star formation at the border of Sh2-294 may 
be due to expansion of HII region. It is therefore necessary to invoke the classical triggering 
scenario of Elmegreen & Lada (1977) that UV radiation compressed the abundant nearby 
gas, leading to the star formation. 



6.2. Age of the ionizing source and the HII region 

The spectral type of ionizing star is of BO ± 0.5 V (see §4.3), which agrees well with 
the spectral type derived from radio continuum emission at 1280 MHz and also from NIR 
CC diagram. Therefore, we assume the spectral type of the massive star of BO V for further 
discussion. The probable age (~ 4x10^ yr) of ionizing star compared to the MS life time 
of BO type star (~ 1 xlO^ yr), implies that the HII region is still under expansion. The 
dynamical age of Sh2-294 can be estimated by using the model by Dyson & Williams (1997), 
assuming ionized gas of pure hydrogen at constant temperature in an uniform medium of 
constant density. If the HII region associated to Sh2-294 is excited by a BO V star, which 
emits ~ 78.75 xlO^^ {LogNiyc ~ 47.90) ionizing photons per second and is expanding into a 
homogeneous medium of density 10^ (10^) cm~^, will form a Stromgren radius (r^) of ~ 0.17 
(0.79) pc. The initial phase is overpressured compared to the neutral gas and will expand 
into the medium of radius r^ at time t as 



n = r.il + ^)'", (4) 

The observed diameter of the ionized gas from GMRT radio contour is ~ 4.7', which 
corresponds to a physical diameter of about 6.6 pc at a distance of 4.8 kpc. An assumed 
sound speed (C//) of 11.4 (4.0) km s^^ in the ionized gas of density 10'^ cm^"^ yields dynamical 
age of ~ 1.5 (4.2) xlO^ yr, whereas in the lower density medium (for example 10^ cm~^) 
the above velocity range will yield ~ 0.4 (1.2) xlO^ yr. However the uncertainty in age 
estimation can not be ignored because of the assumption of a uniform medium and also the 
velocity of expansion changes with time. One would then expect an expanding IF to be 
preceded by a swept-up shell of cool neutral gas as it erodes into a neutral cloud. We have 
detected a half ring of MSX dust emission towards the eastern part of HII region which is 
an indication that neutral gas surrounds the ionized gas. 
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6.3. Triggered steir formation 

Triggered star formation is a sequential process, in which star formation at the interface 
between HII regions and molecular clouds can be triggered by the ionization or shock fronts 
produced by massive stars. The mechanisms mainly responsible for triggered star formation 
have been discussed such as: the "radiation driven implosion" (RDI) and the "coUect-and- 
coUapse" model. In the RDI model, the expanding ionization front compresses the pre- 
existing molecular clumps, leading to density enhancements, which when it exceeds the local 
critical mass, collapses to form new stars (Lefioch & Lazareff 1994; Gorti & HoUenbach 2002). 
In the "collect-and-coUapse" model, the neutral material accumulates between the ionization 
front of the HII region and the shock front in the neutral gas. This compressed, shocked 
material may be dynamically unstable on a short internal-crossing time scale, leading to 
the formation of cometary globules and bright rims (Garcia-segura & Franco 1996) or the 
compressed layer may become gravitationally unstable along its surface, on a long time scale 
(Elmegreen & Lada 1977; Zavagno et al. 2006), leading to the formation of massive fragments 
forming massive stars and/or clusters. 

The CM diagrams (Fig. 8) indicate that ionizing source may have a maximum age of 
~ 4 X 10^ yr. It is rather difficult to estimate the age of the cluster associated with 'region 
A' as it has no optical counterpart and has poor statistics, however the age to the 'region 
A' can be assigned on the basis of the age of CTT stars. The life time of CTTs ranges from 
a few xlO^ to few xlO^ yr. The 'region A' contains two CTT stars and their positions in 
intrinsic (J — H) vs. J CM diagram suggest an age of < 1 x 10^ yr. However, this age 
should be considered as an approximate estimation because of uncertainty in determination 
of extinction as well as error in colors. On the basis of the fraction of NIR excess stars 
(cf. §4.6.1), an age of ~ 1x10^ yr can also be assigned to the 'region A'. If we assume that 
the distance between the ionizing star and the young stars in the 'region A' is closed to a 
projected distance of ~ 1.8 pc (1.3') and that the ionization front expands at a sound speed 
of 4.0 (11.4) km s~^, it would take approximately 4.6 (1.7) x 10^ yr to reach the surface 
of the cloud and can induce second generation stars at the border. The protostellar phase 
in the clouds last for ~ 10^ yr. Therefore, the age difference between the young sources in 
'region A' and central ionizing source, keeping in mind the uncertainty in assigning the ages 
to the clusters, seems to be consistent with the hypothesis that star formation at the border 
may be due to RDI process. 

However, the possibility of "coUect-and-coUapse" process can not also be ruled out by 
considering the following evidence. The radio observation reveals that Sh2-294 is a roughly 
spherically ionized region and is surrounded by half ring of MSX dust emission in the mid-IR 
(see Fig. 17). The ring is enhanced in the middle of the arc of dust emission. We found 
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a MSX point source at the center of dust peak. The presence of dust ring beyond the IF 
indicates that the neutral gas surrounds the HII region. An embedded cluster ('region A') 
is associated with the dust ring within 2MASS detection limit and its central density lies 
close to the position of MSX point source, satisfying most of the predictions proposed by 
the Deharveng et al. (2005) for star formation to occur at the border due to "collect-and- 
coUapse" model. For example the morphology of Sh2-294 region is almost similar to the 
selected HII candidates by Deharveng et al. (2005), which are likely to be examples of 
star formation at the border by this process. However to look for the distribution of dense 
fragments around the Sh2-294, as seen in case of RCW 79 (Zavagno et al. 2006), one needs 
millimeter observations. The fact that the dynamical age of HII region is of ~ 1.5x10^ yr, 
comparable to the age of the stars associated to the 'region A' and the presence of a ZAMS 
B1.5 star further supports the formation of massive star following the collapse of thin shell 
of neutral matter, as predicted by "coUect-and-coUapse" model. 

With the existing observations and data, it will be early to establish or rule out either 
of the scenarios. To provide a conclusive answer to the above mentioned star formation 
scenarios in Sh2-294 we need deep optical and NIR observations, CO, as well as millimeter 
observations. 



7. Summary 

In this paper we provide the first global view of star formation around Sh2-294 region. 
Our multiwavelength data provides a detailed insight into distribution and nature of the 
young stellar objects as well as into the morphology of the thermal ionized gas and dust 
emissions. Based on the 2MASS data, we identify a star cluster at the center of Sh2-294 and 
a distinct group of reddened stars (small cluster) located at 1.8 pc to the east of Sh2-294. The 
optical and near-infrared color-color and color-magnitude diagrams have been constructed 
to estimate the ratio of total-to-selective extinction, reddening, distance and spectral types 
of young stellar objects. The narrow band imaging in Ha and radio continuum map based 
on GMRT observation at 1280 MHz have been used to trace the ionized gas, which show 
interesting morphological details, including an arc-shaped structure in southeast direction 
and creation of lower density zones in the northwest direction through which the ionized gas 
can escape. We identified the most luminous member of the cluster situated at the center of 
the ionized gas and derived its spectral type which is consistent with an ~ BO main sequence 
star. The number of Lyman continuum photons from the radio flux suggests that this single 
star is responsible for the ionization of the nebula. The PDR at the interface between the 
ionized gas and the molecular cloud is traced to the southeast direction by half ring of PAH 
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feature in MSX and diffuse H2 emission. The morphological correlation between ionized gas, 
PDR and H2 emission leads us to conclude that there is a high degree of interaction of UV 
photons from BO star with the nearby molecular cloud in north eastern direction. The small 
cluster is found at the peak of dust ring. The spatial distribution of the YSOs (Class 11 
sources, sources detected in H and K^, bands and only in band) indicates that the region 
around a small cluster is less evolved and suggests that the star formation activity observed 
at the border of Sh2-294 region is possibly triggered by the expansion of HII region. 
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Table 1: Log of observational data. 



Date (UT) 


Filter 


No of 


Exposure (sec) 






frames 


per frame 


2005 Dec 26 


U 


3 


1200 


2006 Mar 06 


U 


2 


300 


2005 Dec 26 


B 


3 


600 


2006 Mar 06 


B 


2 


240 


2005 Dec 26 


V 


3 


600 


2006 Mar 06 


V 


2 


180 


2005 Dec 26 


R 


3 


300 


2006 Mar 06 


R 


2 


120 


2005 Dec 26 


I 


3 


300 


2006 Mar 06 


I 


2 


80 


2006 Sep 28 


Ha 


1 


900 


2006 Sep 28 


Ha cont 


1 


900 


2005 Nov 28 


H2 


21 


90 


2005 Nov 28 




21 


30 



Table 2: Flux density details of the MSX and IRAS point sources 





Flux Deiisit}' ( Jy) for A (//in) 


Source 


/i?A^-HIRES image" 
IRAS PSC 


MSX PSC 




100 60 25 12 


21.3 14.7 12.1 8.3 


IRAS 07141-0920 

MSX6C G224.1880+01.2407 


240.30 125.50 18.03 10.00 
262.40 77.55 7.17 2.82 


3.40 0.86 1.90 0.88 



" Fluxes obtained by integrating over a circular region of diameter 3' centered on the peak. 
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Table 3: Slopes of star distribution in Sh2-294 

Color ratio nT-cluster '^normal 

0.62± 0.02 -0.55 
1.30± 0.04 -1.10 
2.67± 0.07 -1.96 
S.lOi 0.08 -2.42 
3.40± 0.11 -2.60 
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Fig. 1. — (Left) Ha+continuum image (~ 9.7'x 8.6') of the ionized gas associated with 
Sh2-294 region. The box in the Ha image shows the FOV of the NIR image shown in the 
right paneL (Right) Corresponding JHKg color composite image (~ 5.1'x 4.7') from 2MASS 
data. North is up and east is to the left. The circle in the NIR image indicates the young 
embedded cluster at the border of HII region (see text for details). The plus sign in both 
the images marks the position of possible exciting source of Sh2-294 region (see §4.3). 
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Fig. 2. — Surface number density profile as a function of radius for the sources associated 
with Sh2-294 region. The fitted King's profile (dotted line) and the radial variation of 
projected stellar density (solid line) indicates that the cluster is confined within 2.3' radius. 
The error bar represents ±y/N errors. The horizontal dashed line represents the density of 
field stars. 




Fig. 3. — The contour map of the stellar surface number density obtained by counting stars 
in a 12" x 12" (~ 0.28 pc x 0.28 pc) grid for the cluster associated with Sh2-294 region. The 
contours are drawn from 6 to 22 stars/pc^ in steps of 2 stars/pc^. The lowest contour is 
at 3 times of the background level. The positions of the MSX and IRAS point sources are 
marked with asterisk and cross symbol, respectively. The plus sign represents the possible 
exciting star of the region. 
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Fig. 4. — V/V — I CM diagrams for (a) stars in the cluster region (b) stars in the nearby 
field region of same area as of cluster region and (c) statistically cleaned CM diagram, where 
the points with triangles are the probable field stars (see §4.4). 




Fig. 5. — The blue part of the spectrum for the ionizing source of Sh2-294 region. The 
spectrum displays the characteristic similar to BO ± 0.5 V MS member. 
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Fig. 6. — Color-Color {U — B / B — V) diagram for the stars lying within cluster radius. The 
dashed and continuous curves represent intrinsic MS for Z=0.02 by Schmidt-Kaler (1982) 
shifted along the reddening vector of slope 0.72 (arrow line) for E{B — V)mean = 0.2 mag and 
E{B — V)min — 1-35 mag, respectively. The points with triangles are likely to be probable 
field members identified in statistically cleaned CM diagram (see Fig. 4c). 
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Fig. 7. — {V — I) , {V — J) , {V — H) , {V — K) vs. {B — V) two color diagrams for the probable 
cluster members found in statistically cleaned CM diagram. Straight line shows the weighted 
least square fit to the data. 
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Fig. 8. — Field star decontaminated V, B-V and V, V-R CM diagrams for the stars within 
the cluster region. The probable field members identified in Fig. 4c and 6 are not included in 
these figures. The ZAMS in V, B-V and V, V-R CM diagrams (solid curve) by Schmidt-Kaler 
(1982) and Walker (1985) respectively, the isochrone of 4 Myr (dashed curve) by Bertelli et 
al. (1994) and the PMS isochrones for 1 & 5 Myr (dotted curve) by Siess, Dufour, & Forestini 
(2000) are plotted for E{B — V)min = 1-35 mag and distance modulus of 17.6 mag. The open 
circles with dots inside represent the optical counterpart of NIR excess stars (see §4.6.1). 
The arrow points to the possible ionizing source of Sh2-294 region. The average photometric 
errors are shown in the upper right corner of the diagrams. 
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Fig. 9. — {H — K) vs. (</ — H) color-color diagram (a) for the cluster region in Sh2-294 and 
(b) for the control region. The locus of the main sequence (thin solid curve) and giant branch 
(thick dashed curve) are from Bessell & Brett (1988). The dotted line represents locus of 
classical T Tauri stars from Meyer, Calvet, & Hillenbrand (1997). The parallel dashed lines 
are the reddening vectors, with crosses along the lines corresponding to visual extinction 
of 5 mag. The thick continuous curve represents the locus of Herbig Ae/Be stars (Lada & 
Adams 1992). The open squares with dots represent the sources located in 'region A'. The 
plot is classified into three regions namely '•p'-^'-'X" and "P" for different classes of sources 
(see text for details). The average photometric errors are shown in the upper right corner of 
the diagrams. 
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Fig. 10. — Frequency distribution of H — K colors of the sources in the control field (shaded 
area) and the Sh2-294 star forming region (thick solid line). 
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H-K 

Fig. 11. — Color-magnitude diagram {{H — K) vs. K ) for sources detected in all the three 
2MASS bands along with those detected only in H and Ks bands within the cluster radius. 
The thick vertical sohd hues represent the ZAMS curves reddened by Ay=0, 10, 20 mag, 
respectively. The slanting lines trace the reddening vectors for each spectral type. Star 
symbols represent sources of T Tauri type (Class II), filled triangles represent protostellar 
(Class I) objects, crosses represent the probable Herbig Ae/Be stars, filled circles represent 
those sources detected in H and Ks band only and the squares represent the sources which 
lie in 'region A'. The average photometric errors are shown in the upper right corner of the 
diagram. 
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Fig. 12.— (J - H) vs. J intrinsic CM diagram for the YSOs (Class II and I) in Sh2-294 
region. Stars are dereddened according to their individual extinction (see text for details). 
The thin solid curve represents PMS isochrones of 1x10^ yr by Siess, Dufour, & Forestini 
(2000). The arrows represent the reddening vectors for 0.5 and 3.5 Mq PMS stars, respec- 
tively. The symbols are same as in Figure 11. The average photometric errors of the YSOs 
are shown in the upper right corner of the diagram. 
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Fig. 13. — GMRT low resolution map of S1l2-294 region at 1280 MHz. The resolution is 
33" X 17" along PA = 33.3 deg. The contour levels are drawn above three times of rms noise 
of the map, which is 0.26 mJy/beam. 
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Fig. 14. — 2MASS f^s-band image in logarithmic scale overlaid by the GMRT radio contin- 
uum contours at 1280 MHz (blue contours). The plus sign marks the position of the exciting 
source of Sh2-294 region. The 'region A' is shown in the image. 
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Fig. 15. — Comparison of tlie SED from observations and the radiative transfer models of 
the interstellar cloud near 'region A'. The open circles denote IRAS — HIRES data, while 
the crosses represent the MSX data. The HIRES fluxes are integrated over 3' diameter 
circle centred on IRAS peak. The solid curve denotes our best-fit model to the data and 
dashed curve represents the radio- fit model (see the text) . 
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Fig. 16. — The Ha continuum-subtracted contours overlaid over the continuum-subtracted 
H2 (1-0) SI hne image. The prominent peaks in the contours are due to the residuals of 
continuum subtraction around some stars (see Fig. 17). The symobls are same as in Fig. 
14. 
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Fig. 17. — Contour diagram from the MSX A-band image (8.28 /xm), overlaid on an Ha line 
image. The contours start at 2.1x10^ W m~^ Sr~^ and the innermost contour corresponds 
to a value of 2.7x10^ W m^^ Sr^^ Residuals of the continuum subtraction around some of 
the stars are seen in the image. The plus sign represents the position of the central ionizing 
star that emits UV photons. 



Fig. 18. — HIRES processed IRAS maps for Sh2-294 region in the four bands (clockwise 
from the top left)- Tb|^ left) 12 /^m with peak = 41 MJy/Sr, Top right) 25 fxm with peak = 
80.3 MJy/Sr, Bottom left) 60 fim with peak = 384.2 MJy/Sr, and Bottom right) 100 pim 
with peak = 600.5 MJy/Sr. The isophot contour levels in 12 & 25 pim are 95, 90, 80, 70, 
60, 50, 40, 30, 20, 10 & 5%, and in 60 & 100 /im are 95, 90, 80, 70, 60, 50, 40, 30, 20, 10, 5 
& 2.5 % of the respective peaks. The abscissa and the ordinates are in J2000.0 epoch. The 
symbols are same as in Fig. 3. 
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Fig. 19. — Left : The dust optical depth (T25) distribution from the HIRES 12 and 25 
/im maps for the region around Sh2-294 assuming a dust emissivity law of ex oc A~^. The 
contours represent 95, 90, 80, 70, 60, 50, 40, 30, 20 & 10% of the global peak value of 
7.9x10-^. Right : The dust optical depth (tiqo) distribution from the HIRES 60 and 100 
//m maps for the region around Sh2-294 assuming a dust emissivity law of e\ oc A~^. The 
contours represent 95, 90, 80, 70, 60, 50, 40, 30, 20 & 10% of the global peak value of 
1.2x10"^. The symbols are same as in Fig. 3. 
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Fig. 20. — The dust color temperature (T(60/100)) distribution. The contours correspond 
to the temperatures from 30 to 40 K in steps of 2 K. The abscissa and the ordinates are in 
J2000.0 epoch. The symbols are same as in Fig. 3. 
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Fig. 21. — The continuum-subtracted molecular hydrogen 2.12 /im line image from HCT 
overlaid by the GMRT radio continuum contours (blue color). The red contours represent 
the distribution of warm dust from MSX A-band (8.28 /zm). Spatial distributions of Class 
II (star symbols), Class I (filled triangles), probable Herbig Ae/Be sources (crosses), sources 
detected only in H and K^^ bands (filled magenta circles) and sources detected only in 
band (filled red circles) are shown in the figure. 



